Inverse analysis of stochastic moment equations for transient flow in randomly heterogeneous media     

Monica Riva  Alberto Guadagnini  Shlomo P. Neuman  Emanuela Bianchi Janetti  Bwalya Malama 《Advances in water resources》2009

We present a nonlinear stochastic inverse algorithm that allows conditioning estimates of transient hydraulic heads, fluxes and their associated uncertainty on information about hydraulic conductivity (K) and hydraulic head (h  ) data collected in a randomly heterogeneous confined aquifer. Our algorithm is based on Laplace-transformed recursive finite-element approximations of exact nonlocal first and second conditional stochastic moment equations of transient flow. It makes it possible to estimate jointly spatial variations in natural log-conductivity (Y=lnK)$(Y=\ln K)$, the parameters of its underlying variogram, and the variance–covariance of these estimates. Log-conductivity is parameterized geostatistically based on measured values at discrete locations and unknown values at discrete “pilot points”. Whereas prior values of Y at pilot point are obtained by generalized kriging, posterior estimates at pilot points are obtained through a maximum likelihood fit of computed and measured transient heads. These posterior estimates are then projected onto the computational grid by kriging. Optionally, the maximum likelihood function may include a regularization term reflecting prior information about Y. The relative weight assigned to this term is evaluated separately from other model parameters to avoid bias and instability. We illustrate and explore our algorithm by means of a synthetic example involving a pumping well. We find that whereas Y and h can be reproduced quite well with parameters estimated on the basis of zero-order mean flow equations, all model quality criteria identify the second-order results as being superior to zero-order results. Identifying the weight of the regularization term and variogram parameters can be done with much lesser ambiguity based on second- than on zero-order results. A second-order model is required to compute predictive error variances of hydraulic head (and flux) a posteriori. Conditioning the inversion jointly on conductivity and hydraulic head data results in lesser predictive uncertainty than conditioning on conductivity or head data alone.  相似文献   

“Advection,Dispersion and Confusion,”          下载免费PDF全文

Shlomo P. Neuman 《Ground water》2016,54(1):2-3

  相似文献   

Climatology of dust days in the Central Plateau of Iran     

Mesbahzadeh  Tayyebeh  Salajeghe  Ali  Sardoo  Farshad Soleimani  Zehtabian  Gholamreza  Ranjbar  Abbas  Krakauer  Nir Y.  Miglietta  Mario Marcello  Mirakbari  Maryam 《Natural Hazards》2020,104(2):1801-1817

Natural Hazards - Dust storms are a major natural hazard to human health. Severe erosive storms in parts of the Central Plateau of Iran have made the situation very difficult for the inhabitants,...  相似文献   

Editorial: Tribute to Eugene S. Simpson     

Shlomo P. Neuman  Michael E. Campana 《Hydrogeology Journal》1998,6(1):1-2

  相似文献   

Numerical investigation of apparent multifractality of samples from processes subordinated to truncated fBm     

Alberto Guadagnini  Shlomo P. Neuman  Monica Riva 《水文研究》2012,26(19):2894-2908

We investigate numerically apparent multi‐fractal behavior of samples from synthetically generated processes subordinated to truncated fractional Brownian motion (tfBm) on finite domains. We are motivated by the recognition that many earth and environmental (including hydrological) variables appear to be self‐affine (monofractal) or multifractal with Gaussian or heavy‐tailed distributions. The literature considers self‐affine and multifractal types of scaling to be fundamentally different, the first arising from additive and the second from multiplicative random fields or processes. It has been demonstrated theoretically by one of us that square or absolute increments of samples from Gaussian/Lévy processes subordinated to tfBm exhibit apparent/spurious multifractality at intermediate ranges of separation lags, with breakdown in power‐law scaling at small and large lags as is commonly exhibited by real data. A preliminary numerical demonstration of apparent multifractality by the same author was limited to Gaussian fields having nearest neighbor autocorrelations and led to rather noisy results. Here, we adopt a new generation scheme that allows us to investigate apparent multifractal behaviors of samples taken from a broad range of processes including Gaussian with and without symmetric Lévy and log‐normal (as well as potentially other) subordinators. Our results shed new light on the nature of apparent multifractality, which has wide implications vis‐a‐vis the scaling of many hydrological as well as other earth and environmental variables. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Apparent multifractality and scale‐dependent distribution of data sampled from self‐affine processes     

Shlomo P. Neuman 《水文研究》2011,25(11):1837-1840

It has been previously demonstrated theoretically and numerically by the author that square or absolute increments of data sampled from fractional Brownian/Lévy motion (fBm/fLm), or of incremental data sampled from fractional Gaussian/Lévy noise (fGn/fLn), exhibit apparent/spurious multifractality. Here, we generalize these previous development in a way that (a) rigorously subordinates (truncated) fLn to fGn or, in a statistically equivalent manner, (truncated) fLm to fBm; (b) extends the analysis to a wider class of subordinated self‐affine processes; (c) provides a simple way to generate such processes and (d) explains why the distribution of corresponding increments tends to evolve from heavy tailed at small lags (separation distances or scales) to Gaussian at larger lags. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

The maximal runaway temperature of Earth-like planets     

Nir J. Shaviv  Giora Shaviv  Rainer Wehrse 《Icarus》2011,216(2):403-414

In Simpson’s (Simpson, G.C. [1927]. Mem. R. Meteorol. Soc. II (16), 69–95) classical derivation of the temperature of the Earth in the semi-gray model, the surface temperature diverges as the fourth root of the thermal radiation’s optical depth. No resolution to this apparent paradox was yet obtained under the strict semi-gray approximation. Using this approximation and a simplified approach, we study the saturation of the runaway greenhouse effect.First we generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength domain of the absorption.Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the greenhouse effect saturates and that any Earth-like planet has a maximal temperature which depends on the type of and distance to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the visible, could have at most a surface temperature of 1200–1300 K irrespective of the nature of the greenhouse gas.We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the value of λ_cut, the wavelength above which the atmosphere becomes optically thick. Unless λ_cut is small and resides in the optical region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. This has an interesting observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, saturation takes place by radiating through windows in the FIR.  相似文献   

Non-local and localized analyses of non-reactive solute transport in bounded randomly heterogeneous porous media: Theoretical framework     

Eric Morales-Casique  Shlomo P. Neuman  Alberto Guadagnini 《Advances in water resources》2006

  相似文献   

Longitudinal dispersivity data and implications for scaling behavior     

Neuman SP 《Ground water》2006,44(2):139-140

  相似文献   

The theory of steady-state super-Eddington winds and its application to novae     

Nir J. Shaviv 《Monthly notices of the Royal Astronomical Society》2001,326(1):126-146

We present a model for steady-state winds of systems with super-Eddington luminosities. These radiatively driven winds are expected to be optically thick and clumpy as they arise from an instability-driven porous atmosphere. The model is then applied to derive the mass loss observed in bright classical novae. The main results are as follows.
(i) A general relation between the mass-loss rate and the total luminosity in super-Eddington systems.
(ii) A natural explanation of the long-duration super-Eddington outflows that are clearly observed in at least two cases (Novae LMC 1988 #1 and FH Serpentis).
(iii) A qualitative agreement between the prediction and observations of the mass loss and temperature evolution.
(iv) An agreement between the predicted average integrated mass loss of novae as a function of white dwarf mass and its observations.
(v) A natural explanation for the 'transition phase' of novae.
(vi) Agreement with η Carinae, which was used to double check the theory: the prediction for the mass shed in the star's great eruption agrees with observations to within the measurement error.  相似文献